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Abstract. It is shown that any singular Lagrangian theory: 1) can be for- 
mulated without the use of constraints by introducing a Clairaut-type version 
of the Hamiltonian formalism; 2) leads to a special kind of nonabelian gauge 
theory; 3) coincides with the many-time classical dynamics. A generalization 
of the Legendre transform to the case when the Hessian is zero is done by 
using the mixed (envelope/general) solution of the multidimensional Clairaut 
equation. The corresponding system of equations of motion is equivalent to 
the Lagrange equations and has a linear algebraic subsystem for "unresolved" 
velocities. Then the equations of motion are written in the Hamiltonian form 
by introducing a new bracket. This is a "shortened" formalism, since initially 
it does not contain "nondynamical" (degenerate) momenta at all, and there- 
fore there is no notion of constraint. It is shown that any classical degenerate 
Lagrangian theory in its Clairaut-type Hamiltonian form is equivalent to the 
many-time classical dynamics. 
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1. Introduction 

Nowadays, many fundamental physical models are based on gauge field theo- 
ries [T31ISS]. On the classical level, they are described by singular (degenerate) 
Lagrangians, which makes the passage to a Hamiltonian description, which is im- 
portant for quantization, highly nontrivial and complicated (see, e.g., |45 p 51j). 

A common way to deal with such theories is the Dirac approach |15) based on 
extending the phase space and constraints. This treatment of constrained theories 
has been deeply reviewed, e.g., in lecture notes [S^] and books [23,t30j. In spite 
of its general success, the Dirac approach has some problems [16l|43l[57] and is 
not directly applicable in some cases, e.g., for irregular constrained systems (with 
linearly dependent constraints) ^6f,41 . 

Therefore, it is of considerable interest to reconsider basic ideas of the Hamilton- 
ian formalism in general from another point of view. The approach presented here 
does not depend on the structure of constraints and hopefully will be applicable to 
all cases. 

In the standard approach for nonsingular theories [4j|49], the transition from 
Lagrangian to Hamiltonian description is done by the Legendre transform, and 
then finding the Hamiltonian as an envelope solution of the corresponding Clairaut 
equation |32| . The main idea of our formulation is the following [20': for singular 
theories, instead of the Lagrange multiplier procedure developed by Dirac jl5j . 
alternatively we construct and solve the corresponding multidimensional Clairaut 
equation |32j . In this way, we state that the ordinary duality can be generalized to 
the Clairaut duality. 

In this paper we develop our previous work [T51[^ to construct a self-consistent 
analog of the canonical (Hamiltonian) formalism and present a general algorithm 
to describe any Lagrangian system (singular or not) as a set of first-order differ- 
ential equations without introducing Lagrange multipliers. The connection of our 
approach with the Dirac constraints method is shown in Appendix |B] 

The advantage of this "shortened" formalism, when the "nondynamical" (de- 
generate or "non-physical") momenta do not enter in initial consideration, is its 
clarity and simplicity. Similar formulations (without primary constraints) were 
treated in p4!f25{l38], while here we explain the reason and source of constraints 
appearance from another viewpoint. 

The "shortened" approach can play an important role for quantization of such 
complicated constrained systems as gauge field theories _34] and gravity ^44] . These 
questions will be considered in future papers. To illustrate the power and simplicity 
of our method, we consider such examples, as the Cawley Lagrangian [91 (which 
led to difficulties in the Dirac approach), and the relativistic particle. A novel 
Hamiltonian-like form of the equations of motion is achieved by defining a new 
(non-Lie) bracket which is not anticommutative and does not satisfy the Jacobi 
identity. Note that quantization of such brackets can be done by means of non-Lie 
algebra methods (see, e.g., [42l[52]). This will be a subject of a subsequent paper. 

While analyzing the part of the equations of motion corresponding to "unre- 
solved" velocities, we arrive effectively at a kind of nonabelian gauge theory in the 
"degenerate" coordinate subspace, related to these velocities. Here we have shown 
that the Clairaut-type formulation is equivalent to the many-time classical dynam- 
ics developed in [I7l[36] if "nondynamical" (degenerate) coordinates are treated 
as additional "times" . Finally, in an Appendix we show that, after introduction 
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of "nondynamical" momenta, corresponding Lagrange multipliers and respective 
constraints, the Clairaut-type formulation presented here corresponds to the Dirac 
approach [IP. 

In general, the paper has also a clear methodological aspect: it shows that one 
can deal with well-known problems from a different than commonly used viewpoint. 
This can lead to better understanding and new insights into such an important 
subject as Hamiltonian formalism in itself. 

To simplify matters, we use coordinates, but all the statements can be readily 
reformulated in a coordinate free setting [5J[n3]. We consider systems with a finite 
number of degrees of freedom. The Clairaut-type formulation of field theories will 
be done elsewhere. 

2. The Legendre-Fenchel and Legendre transforms 

We start with a brief description of the standard Legendre-Fenchel and Legendre 
transforms for the theory with nondegenerate Lagrangian [3147]. Let0L(g^,w^), 
A = 1, . . .n, be a Lagrangian given by a function of 2n variables (n generalized 
coordinates and n velocities — = dq-^/dt) on the configuration space 
TM, where M is a smooth manifold. We consider the time-independent case for 
simplicity and conciseness, which will not infiuence the main procedure. 

By the convex approach definition (see e.g. [5ll46] ). a Hamiltonian H is a 

dual function on the phase space T*M (or convex conjugate [47]) to the Lagrangian 
(in the second set of variables pa) and is constructed by means of the Legendre- 

fi Fen 

Fenchel transform L <^ H^"'' defined by [HllS] 

H^°^ {q^,PA) ^ sup G{q^,v^, pa), (2.1) 

yA 

n 

G {q^, v^,pa) = E PB^"" ~ L (<Z^, v^) . (2.2) 

B=l 

Note that this definition is very general, and it can be applied to nonconvex 
[2] and nondifferentiable [56] functions L(^q'^,v^), which can lead to numerous 
extended versions of the Hamiltonian formalism (see, e.g., [12l[3Tl|48]). Also, a 
generalization of the convex conjugacy can be achieved by substituting in (|2.2p the 
toTm pav^ by any function vjf (^pA,v"^) satisfying special conditions 26 . 

In the standard mechanics [29], one usually restricts to convex, smooth and 
differentiable Lagrangians (see, e.g., [5]|50] ). Then the coordinates q'^{t) are treated 
as fixed (passive with respect to the Legendre transform) parameters, and the 
velocities v^{t) are assumed independent functions of time. 

According to our assumptions the supremum (|2.1|) is attained by finding an 



extremum point = w^j^ of the ( "pre-Hamiltonian" ) function G (q"^, u"^,p^) 



extr 

which leads to the supremum condition 

dL (q^,v^) 

PB 



(2.3) 



It is commonly assumed (see, e.g., [Sl[5ni|^) that the only way to get rid of 
dependence on the velocities in the r.h.s. of (|2.ip is to resolve (|2.3p with respect 



^We use indices in arguments, because we need to distinguish different kinds of coordinates 
(similar to [55)). For the same reason, we use the summation signs with explicit ranges. 
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to velocities and find its solution given by a set of functions 

vg,r^V^{q^,PA). (2.4) 

This can be done only in the class of nondegenerate Lagrangians L(g'^,w'^) — 
^nondcg .yA^j ^-j^^ -j-j^g sccond sct of Variables v^), which is equivalent to the 
Hessian being non-zero 



det 



g2^nondeg (g^^^;^) 



^0. (2.5) 



Then substitute w^j^ to (|2.1I) and obtain the standard Hamiltonian (see, e.g., [5i29) ) 

H {q^,PA) ''^ G {q^, vttr,PA) 

71 

The passage from the nondegenerate Lagrangian i"°ndcg ^^a^^a^^ ^^i^ Hamil- 
tonian H(^q^,pA) is called the Legendre transform which will be denoted by 

^nondcg ~'S^ 

From the geometric viewpoint [Tl l39[[51] . the Legendre transform being applied 
to L (^q^,v^) is tantamount to the Legendre transformation from the configuration 
space to the phase space Leg : TAf T*M (or between submanifolds in the 
presence of constraints 



3. The Legendre-Clairaut transform 

An alternative way to deal with the supremum condition (j2.3p is to consider 
the related multidimensional Clairaut equation [19j . The connection between the 
Legendre transform, convexity and the Clairaut equation has a long story [331149] 
(see also [1]). 

We differentiate (|2.6p by the momenta pA and use the supremum condition 
to get 



dH{q^,PA) 




PC 



^ {q\pA) 
dL {q^,v^ 



dV^iq^^PA, ^,s,A 



dpB 



V''{q^,PA), (3.1) 



which can be called the dual supremum condition (indeed this gives the first set of 
the Hamilton equations, see below). The relations (|2.3p . (j2.6p and (|3.ip together 
represent a particular case of the Donkin theorem (see e.g. [29j). 
Then we substitute (|3.1I) in (12. 6p and obtain 



"'^ 1 1 A I . f , t / I I I 1 1 . I I A I \ 

OpB \ Opc 

which contains no manifest dependence on velocities at all. It is important that, 
for nonsingular Lagrangians, the relation p.2p is an identity, which follows from 
and p.ip by our construction. 
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Now we make the main step: to treat the equation p.2p by itself (without refer- 
ring to (|2.3|) . (|2.6p and p.ip ) as a definition of a new transform being a solution of 
the foUowing nonhnear partial differential equation (the multidimensional Clairaut 
equation) [19^ 



B=l 



dpB 



-Ll q 



dpA 



(3.3) 



where pA are parameters (initially not connected with pA defined by (|2.3p ) and 
L ((j"^, va) is any differentiable smooth function of 2n variables, and here we do not 
demand the nondegeneracy condition (|2.5p . We call the transform defined by 

-2<^9, TrCl 



L H 



and 



a Clairaut duality transform (or the Leg endre- Clairaut transform) 
Pa) a Hamilton- Clairaut function |19j . 
Note that (|2.3|) is commonly treated as a definition of all dynamical momenta 
PA, but we should distinguish them from the parameters of the Clairaut duality 
transform pA- In our approach, before solving the Clairaut equation and apply- 
ing the supremum condition (12. 3p . the parameters pA are not connected with the 
Lagrangian. Thus, the corresponding Leg endre- Clairaut transformation is the map 



Leg^' : TM 

and we call T^M a Clairaut (extended) phase space. 

The difference between the Legendre- Clairaut transform and the Legendre trans- 
form is crucial for degenerate Lagrangian theories |19) . Specifically, the multidi- 
mensional Clairaut equation p.3|) has solutions even for degenerate Lagrangians 



where the space T 



ci 



M in local coordinates is just iq^. 



pa) 



(3.4) 

is another (along with 
in the case of degen- 



L [q^^ v^) = L^°^ (^q-^, w"^) when the Hessian is zero 



det 



= 0. 



In this way, the Legendre-Clairaut transform (13.31) £eg 
Hcfl^''") counterpart to the ordinary Legendre transform 
erate Lagrangians. 

To make this manifest and to find solutions of the Clairaut equation 
differentiate it by pc to obtain 



E 

5=1 



PB 



dL 



d^H^' [q^^PA) 

dpBdpc 



we 



(3.5) 



Now we apply the ordinary method of solving the Clairaut equation (see Appendix 
E)) . There are two possible solutions of p.5l) . one in which the square brackets van- 
ish (envelope solution) and one in which the double derivative in velocity vanishes 
(general solution). The l.h.s. of p.Sp is a sum over B and it is quite conceivable 
that one may vanish for some B and the other vanish for other B. The physical 
reason of choosing the particular solution is presented in Section HI Thus, we have 
two solutions of the Clairaut equation: 

1) The envelope solution defined by the first multiplier in (13. 5p being zero 



PB = PB = 



dL {q^,v^) 



(3.6) 
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which coincides with the supremum condition ()2.3|) . together with p.ip . In this 
way, we obtain the standard Haniiltonian ()2.6p 

H?L (q^.PA) |p.=p. = H {q^,pA) . (3.7) 
Thus, in the nondegenerate case, the "envelope" Legendre-Clairaut transform 
■^^flenu • ^ ~^ -^ent) coincides with the ordinary Legendre transform by our con- 
struction here. 

2) A general solution is defined by 



dpsdpc 



(3.8) 



which gives ^ — ^Pa) _ ^ here c^ are arbitrary smooth functions of q'' 

OPB 

and the latter are considered in the Clairaut equation p.3p as parameters (passive 
variables). Then the general solution acquires the form 

iJ«„(q^,p^,c^) = f^p5c^-L(<z^,c^), (3.9) 

which corresponds to a "general" Legendre-Clairaut transform 2tQ^l„ : L — > H^^„. 
Note that the general solution i?^g„ (q"^,pA,c"^) is always linear in the variables 
PA and the latter are not actually the dynamical momenta pA, because we do not 
have the envelope solution condition (|3.6p , and therefore now there is no supremum 
condition (|2.3p . The variables are in fact unresolved velocities in the case of 
the general solution. 

Note that in the standard way, S-iQenv ^an be also obtained by finding the enve- 
lope of the general solution , i.e. differentiating p.9[) by as 

-^^ — d^-^ ^^-^"^ 

which coincides with p.6p and \2.2>\ . This means that H^^n {i^^pA, c^) |c^=t,^ is in 
fact the "pre-Hamiltonian" (j2.2p , which was needed to find the supremum in (|2.ip . 
Let us consider the classical example of one-dimensional oscillator. 

Example 3.1. Let L (x,w) = tou^/2 — kx'^ /2 (m, k are constants), then the corre- 
sponding Clairaut equation p.3p for H = iJ^' (a;,p) is 

m , ,s2 kx"^ 
H^pH'^~-{H'^) +— , (3.11) 

where prime denotes partial differentiation with respect to p. The general solution 
is 

H^L i.x,p,c)^pc^ — + —, (3.12) 

where c is an arbitrary function ("unresolved velocity" v). The envelope solution 
(with p — p) can be found from the condition 



do 

which gives 



in the standard way. 



n , P 

p - mc = Cextr = — , 

m 



HZAx,p)^^ + ^ (3.13) 
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Example 3.2. Let L (x, v) — xexp kv, then the corresponding Clairaut equation for 
H = H'^^ is 

H = pH'p-x exp (kH'p) . (3.14) 

The general solution is 

H^en i^: p) = pc — X exp kc, (3.15) 
where c is any smooth function of x. 

The envefope solution (with p — p) can be found by differentiating the general 
solution (I3.15P 



ci 



dc 



p — X exp fee — 



k X 



which leads to 



H?nv {x,p) = |ln^ - p. 
4. The mixed Legendre-Clairaut transform 



Now consider a singular Lagrangian L (^q^ 



(3.16) 



for which the 



Hessian is zero (|3.4p 



This means that the rank of the Hessian matrix Wab — 



is r < n, and we suppose that r is constant. We rearrange indices 
of Wab in such a way that a nonsingular minor of rank r appears in the upper 
left corner [24]. Represent the index A as follows: if A = 1, . . . ,r, we replace A 
with i (the "regular" index), and, if A — r + 1, . . . , n we replace A with a (the 



"degenerate" index). Obviously, detWij ^ 0, and rankWij 



Thus any set of 



variables labelled by a single index splits as a disjoint union of two subsets. We call 
those subsets regular (having Latin indices) and degenerate (having Greek indices) . 

The standard Legendre transform £eg is not applicable in the degenerate case 
because the condition p.Sp is not valid [8l[53]. Therefore the supremum condition 
(|2.3p cannot be resolved with respect to degenerate A, but it can be resolved only 
for regular A, because deiWij ^ 0. On the contrary, the Clairaut duality trans- 
form given by p.3|) is independent of the Hessian being zero or not [121 ■ Thus, 
we state the main idea of the formalism we present here: the ordinary duality can 
he generalized to the Clairaut duality. This can be rephrased by saying that the 
standard Legendre transform £eg (given by (j2.6|) ) can be generalized to the singu- 
lar Lagrangian theory using the Legendre-Clairaut transform £eg^' given by the 
multidimensional Clairaut equation p.3p . 

To find its solutions, we differentiate p.3|) by pA and split the sum (|3.5p in B as 
follows 



E 

4=1 



E 



dL {q^,v^) 
dL (g^. 



dpidpc 



PA 



dpadpc 



As det Wij 7^ 0, we suggest to replace (|4.1|) by the conditions 



Pi = Pi 



dL {q^,v^) 



L , . . . , I , 



d^H^' {q^,PA) 

dpadpc 



0, a = r + 1, . . .n. 



(4.1) 

(4.2) 
(4.3) 
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In this way we obtain a mixed envelope/general solution of the Clairaut equation 
as fohows [19 . After resolving (|4.2p by regular velocities v"^ — y (g"^,pi,c") and 

writing down a solution of (14. 3p as — ^ — c" (where c" are arbitrary 

functions), we obtain a mixed Hamilton- Clairaut function 

r 

n 

+ ^ (4.4) 



which is the desired "mixed" Legendre-Clairaut transform L i— ii^" H^^^ written 
in coordinates. 

Note that (|4.4p coincides with the "slow and careful Legendre transformation" 
of [55] and with the "generalized Legendre transformation" of ^Oj, while we have 
obtained in a new way. 

Example 4.1. Let L [x, y, Vx, Vy) = myv'^/2 + kxVy, then the corresponding Clairaut 
equation for H = i/*-" {x,y,px,Py) is 

H = PxH;^ + pyH^^ -!!^(Hif- kxH'^^ . (4.5) 

The general solution of (|4.5I) is 

2 

rrCl / -- I- _ '"^y^x _ 1 

^ gen Px 1 Py T Cxt Cy ) — PxCx ~r PyCy ^ KXCy, 

where Cx, Cy are arbitrary functions of the passive variables x, y. Then we differen- 
tiate 

nrrCl 

^^qen ^ Fxtr Px 

—2— =px- mycx = 0, =^ c^^"^ = , 

OCx my 

^ gen _ , 

dCy 

Finally, resolve only the first equation and set Cy i — > Vy an "unresolved velocity" 
to obtain the mixed Legendre-Clairaut transform (cf. [55], Examples 5 and 17) 

H^\x ix,y,Px,Py,Vy) = iPy " M ■ (4-6) 

5. Hamiltonian formulation of singular Lagrangian systems 

Let us use the mixed Hamilton-Clairaut function H^^^ (^q^ ,pi,pa,v°') (|4.4p to 
describe a singular Lagrangian theory by a system of ordinary first-order differential 
equations. In our formulation we divide the set of standard Lagrange equations of 
motion 

d dL{q^,v^) _ dL{q^,v^) 

dt dv^ dqB ^ ^ 

into two subsets, according to the index B being regular {B — i — l,...,r) or 
degenerate {B — a ^ r + 1, . . . n). We use the designation of "physical" momenta 



GENERALIZED DUALITY, HAMILTONIAN FORMALISM AND NEW BRACKETS 9 



(|4.2p in the regular subset only, such that the Lagrange equations become 

dpi dL {q^ 



dt 



dBa {q ,Pi 
di 



dL {q^,v' 



where 



def 



dL iq' 



(5.2) 
(5.3) 

(5.4) 

are given functions which determine dynamics of the singular Lagrangian system 
in the "degenerate" sector. The functions Ba {q^,pi) are independent of the unre- 
solved velocities since rankHOiB = r. One should also take into account that 
now 

dt 

Note that before imposing the Lagrange equations (|5.2p (while solving the Clairaut 
equation (IB.Sp ). the arguments of L (q"^, v^) were treated as independent variables. 

A passage to an analog of the Hamiltonian formalism can be done by the standard 
procedure: consider the full differential of both sides of ()4.4|) and use the supremum 
condition ()4.2|) . which gives (note that in previous sections the Lagrange equations 
of motion (I5.ip were not used) 



dt 



(5.5) 



dH. 



ci 



dpt 



V'{q^,p.„v'') 



dH, 



Cl 



dpa 
Cl 



dH, 



dL {q 



dq 

dH, 



Cl 

mix 



dq' 
dL (q^, 



I3=r+l 



Bp{q^.p.:)] 



dq° 



dq°' 



dvl^ 



dvP 



An application of (|5.2p yields the system of equations which gives a Hamiltonian- 
Clairaut description of a singular Lagrangian system 



(5.6) 
(5.7) 
(5.8) 



^ mix 


dq' 


dpi 


~ H' 


dH^y 

mix 


dq" 


dpa 


~ ~dt' 


dH^^. 

^ mix 


_ dpi 


dq' 


~ dt 


dH"^^ 

^ mix 


dB^{ 



f E [PP-Bp{q^,p^)] 
A \ " 



dvi^ 



dv^ 



dt 



P=r+1 



[pp-Bp — 



(5.9) 



The system (|5.6p - (|5.9p has two disadvantages: 1) it contains the "nondynamical" 
momenta Pa] 2) it has derivatives of unresolved velocities w". We observe that 
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we can get rid of these difficulties, if we reformulate (|5.6|) - (|5.9|) by introducing a 
"physical" Hamiltonian 

n 

Hphys {q^,P^) = {q^,P^,Pa,V^) " ^ [p^ " 5/3 {q^,P^)] (5.10) 

/3=r+l 

which does not depend on the degenerate variables ■ The "physical" Hamiltonian 
(|5.10p can be rewritten in the form 



H,Hys {q^,P^)^Y.P^V^ {q^,P^,v'^) 

n 

+ J2 B„{q^,p,)v^-L{q^,V'{q^,p,,v^),v"), (5.11) 

a— r+l 

Then using (14. 2p . one can show that the r.h.s. of (j5.10p indeed does not depend 
on "nondynamical" momenta pa and degenerate velocities v" at all (which justifies 
the term "physical"). Hereafter we will use the shorthand Bp := Bjs {q^,Pi), but 
be sure to remember it depends on (^q^,pij. 
Now we introduce a "(/"-long derivative" 

dX 

where X — X [q^,Pi) is a smooth scalar function, and 

V— \ / dX dY dY dX \ 

is the physical Poisson bracket (in regular variables q*, p, ). 

Then we obtain from (|5.6p - (l5.9p the main result of our Clairaut-type formulation: 
the sought- for system of ordinary first-order differential equations {the Hamilton- 
Clairaut system of equations) which describes any singular Lagrangian classical 
system has the form 

'i-{q\H,,ys{q\p.)],,,^- E {q\B,]^^J-^, ^^l,...r (5.14) 

/3=r+l 



" dq^ 
E Fa0 {q^,Pi) ^ DaHphys {q^,Pi) , a^r + l,...,n. (5.16) 

We also introduce in (|5.16p a "q"-field strength" of the "g"-gauge fields" 

dB 3B 

which is nonabelian due to the presence of the Poisson bracket in the physical phase 
space. 

The system (j5.14p - (l5.16l) is equivalent to the Lagrange equations of motion (15. ip 
by construction. Thus, the Clairaut-type formulation (|5.14p - (|5.17p is valid for any 
Lagrangian theory, as opposite to other approaches. A nonsingular system contains 
no "degenerate" variables at all, because the rank of the Hessian r is full {r — n). 
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The distinguishing property of any singular system (r < n) is clear and simple: it 
contains the additional system of linear algebraic equations ()5.16p for "unresolved" 
velocities which can be analyzed and solved by standard linear algebra methods. 

Indeed the system (I5.16[) gives a full classification of all singular Lagrangian 
theories, which is done in the next section. 

Example 5.1 (Cawley [5]). Let L = xy + zy'^/2, then the equations of motion are 

x = yz, 2/ = 0, y^=0. (5.18) 

Because the Hessian has rank 2 and the velocity z does not enter into the La- 
grangian, the only degenerate velocity is i {a — z), the regular momenta are 
Px = y, Py ^ X (i ^ x, y). Thus, we have 

1 2 

Hphys = PxPy - -j^zy , B^=0. 
The equations of motion (|5.14p - (l5.15l) are 

Px = 0, py= yz, (5.19) 

and the condition (j5.16p gives 

D^Hphys^^^^-ly'^0. (5.20) 

Observe that (|5.19|) and (|5.20[) coincide with the initial Lagrange equations of mo- 
tion ([ET5)) . 

Example 5.2. The Classical relativistic particle is described by 



L = -mR, R= il- if, (5.21) 



i—x^y^z 



where a dot denotes derivative with respect to the proper time. Because the rank 
of the Hessian is 3, we consider velocities ii as regular variables and the velocity 
Xq as a degenerate variable. Then for the regular canonical momenta we have 
Pi — mXi/R which can be resolved with respect to the regular velocities as 



y i—x,y,z 



Using (|5.4p and (|5.11l) we obtain for Hamiltonians 

Hphys = 0, h.j;g = = -E, (5.23) 

H 

and so the "physical sense" of h^gis that it is indeed the energy (|5.22p . Equations 
of motion (I5.14p - (|5.15p are 

. P» . dh^g . 
Xt^Xo — , Pi = ^ — xo^O, 
E oxi 

which coincide with the Lagrange equations following from (|5.2ip directly, the ve- 
locity xq is arbitrary. 
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6. NONABELIAN GAUGE THEORY INTERPRETATION 

Let US consider "q^-long derivative" ()5.12|) and "g^-field strength" (|5.17|) in more 
detail. Note that the "g^-long derivative" satisfies the Leibniz rule 

{Bp,B,}^^y^ = {D^Bp,B,}^^^^^ + {Bp,D^B^}^^y^ (6.1) 

which is valid while acting on "g"-gauge fields" Ba ■ The commutator of the "g"- 
long derivatives" is now equal to the Poisson bracket with the "(/"-field strength" 

{D^Dp-DpD^)X^{F^p{q^,p;),X]^^^^. (6.2) 
In the particular case, acting on Fap {q^,pi), this gives 

D^DpF^p{q^,p,) ^Q. (6.3) 
Let us introduce the i3a-transformation 

^B„X = {i?a,XU,,, (6.4) 

which satisfies 

{^BjBf, ~ Sb^SbJ) B^ = 8{B^,Bf:}^^^^B^. (6.5) 

Sb^F^-, {q^.Pr) = {D^Dp - DpD^) B„, (6.6) 
Sb^ {Bp,B,}^^^^ = {6b^Bp,B,}^^^^ + {B0,6BMp,ys ■ (6-7) 

This means that the "g^-long derivative" Da (|5.12p is in fact a "g"-covariant de- 
rivative" with respect to the Ba-transformation (|6.4|) . Indeed, observe that "Dq 
transforms as fields" ()6.4|) . which proves that it is really covariant (note the cyclic 
permutations) 

^B^DpBj + SB^DaBp + SB^DjBa 

^ {Ba,DpB,}^^^^ + {B„ DM^^^^ + {B0,DM^,^^^ . (6.8) 

The "g"-Maxweir' equations of motion for the "(/"-field strength" are 

DaFap{q'^,Pr)^j0{q^,Pr), (6.9) 
DaFp^ {q'',p,) + D^Fap {q^,p^) + DpF^^ {q^,P^) = 0, (6.10) 

where Ja {q^,Pi} is a "(/"-current" which is a function of the initial Lagrangian ()2.2p 
and its derivatives up to third order. Due to (j6.3p the "g"-current" is conserved 

DaJa{q'^,P^) =0. (6.11) 

Thus, a singular Lagrangian system leads effectively to some special kind of the 
nonabelian gauge theory in the "degenerate" coordinate subspace q", in which 
"nonabelianity" appears not due to a Lie algebra, as in the Yang-Mills theory, 
but "classically", due to the Poisson bracket in the physical phase space {q^,pi). 
The corresponding manifold can perhaps be interpreted locally as the degenerate 
Poisson manifold being a direct product of real space of dimension (n — r) and 
symplectic manifold of dimension (r, r) , where r is the rank of Hessian. 
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7. Classification, gauge freedom and new brackets 
Next we can classify singular Lagrangian theories as follows: 

(1) Gaugeless theory. The rank of the skew-symmetric matrix Fap {q^,Pi) 
is "full", i.e. lankFap {q^,Pi) — n — r and constant, and so the matrix 
Fap {q^,Pi) is invertible, and all the (degenerate) velocities v" can be found 
from the system of linear equations ()5.16p in a purely algebraic way. 

(2) Gauge theory. The skew-symmetric matrix Fai3 {q^,Pi) is singular. If 
rank^Fo,^ [q^,Pi) = rj < n — r, then a singular Lagrangian theory has 
n — r — rf gauge degrees of freedom. We can take them arbitrary, which 
corresponds to the presence of some symmetries in the theory^ Note that 
the rank rf is even due to skew-symmetry of Fap (^q^,pi). 

In the first (gaugeless theory) case one can resolve (|5.16|) as follows 

n 

FP''{q^.P^)D^H^hys{q^,P^), (7.1) 

a— r+1 

where -F"^ {o.^^Vi) is the inverse matrix to -Fq,/:; {Q^^Pi}^ 

n 71 

J2 P^) F^' - E (^^'^'O i'i^^P^) = ^l- (7.2) 

/3=r+l P=r+1 

Substitute (|7.1I) in (|5.14l) - (|5.15p to present the system of equations for a gaugeless 
degenerate Lagrangian theory in the Hamiltonian-likc form as follows 

%-{l\H,,ys{q^P.)},^^^, (7.3) 
^-{p^,H,,ys{q^p.)}^_, (7.4) 
where the bracket is defined by 

n n 

{X,Y}^e^ = {X,Y}^,^^- J2 E {X,B^}^,^J^^{q^,p,)DpY. (7.5) 

a=r+l fS=r+l 

Then the time evolution of any function of dynamical variables X {q^ ,pi) is also 
determined by the bracket (|7.5p as follows 

^ = {X,iI,,,.(g^pO}„e.- (7-6) 

In the second (gauge theory) case, with the singular matrix Fap {q^,Pi) of rank 
r f , we rearrange its rows and columns to obtain a nonsingular rf x rf submatrix in 
the left upper corner. In such a way, the first rf equations of the system of linear 
(under also rearranged v^) equations (j5.16p are independent. Then we express the 
indices a and /3 as pairs a = (ai, a2) and /? = (/3i, /32), where ai and /3i denote the 
first rf rows and columns, while a2 and /32 denote the rest of n — r — rf rows and 



^Here we do not analyze these symmetries, which will be the subject of a separate paper. 
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columns. Correspondingly, we decompose the system (|5.16p as 

r+rf n 
I3i=r+1 h=r+ri + l 

(7.7) 

f)l=r+l I32=r+rf + l 

(7.8) 

Because the matrix Fa^p-^ {q^,Pi) is nonsingular by construction, we can find 
the first r/ velocities 

r+rf 

yP^^ Y F^^'^^ {q^,p,)D^,Hp^,y4q^,p,) 

ai — r+l 

- Y F^'"' (.q^,P^)Fo.,p,{q^,p,)v^^ (7.9) 

ai— r+l 

where F^'-"^ {l^jPi) is the inverse of the nonsingular r/ x r/ submatrix 
^ai/Ji {q^,Pt) satisfying 

Then, since ranki^^^ {q'^jPi) = i"f, the last n — r — rf equations (|7.8I) are linear 
combinations of the first r/ independent ones (j7.7p . which gives 

r+rf 

Fo.,^,{q^,P^)^ Y Kl{q^,P^)Fo.,P,{q^,P^), (7.10) 
ai — r+l 

Fc.hiq^.P)^ Y K\{q^.P^)Fo.,0,{q^.P^), (7.11) 
ai— r+l 
r+rf 

Do.,Hphys{q^,P^^^ Y K\{q^.P^)Do.,Hphys{q^.P^), (7.12) 
Qi— r+l 

where A"^ {q'^ ,Pi) are some r/ x [n — r — rf) smooth functions. Using the relation 
(|7.10p and invertibility of Fq^^^ {q^,Pi) we eliminate the functions (9"^, Pi) and 
obtain 

r+r^' r+Tf 

Fo.,p,{q^.Pr)^ Y E ^O2/^i(9^,k)^^'"'('?'^'P0^"i&('?^'P0' 
Qi=r+1 ^i=r+l 

(7.13) 

Dc.,Hphy,{q^,p,)^ Y Y F^,^,{q^,p,)F^'''' {q^,p,)D^,Hphys{q^,P^)■ 

ai— r+l ^1— r+l 

(7.14) 

This indicates that a gauge theory is fully determined by the first rf rows 
of the (rearranged) matrix Fap {q^,Pi) and the first rf (rearranged) derivatives 
DaiHphys {q^,Pi} only. Next, using further linear combinations of full rows, we 
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can set all elements Fa-^is^ {q^,Pi) — 0, and obtain 

Da^Hphys {q'^^Pi) ^0, a2 = r + r/ + 1, . . . ,n, (7-15) 
which means that functions Da2Hphys {q^TPi) manifestly do not depend on the 

Thus, in addition (to the gaugeless case), now we also have n — r — r f relations 
()7.15p (which are, in fact, combinations of motion integrals for the system (15. 14^ - 
(|5.15p V We can make the unresolved n — r — r/ velocities vanish 

v'^^ = (7.16) 

by some "gauge fixing" condition. Then (|7.9I) becomes 

r+rj 

^ F^^"^{q^,p,)Do,,HpUys{q^,P^)■ (7.17) 
ai— r+1 

By analogy with (|7.3p - (|7.4p . we can also write the system of equations for a singular 
Lagrangian theory in the Hamiltonian form (in the gauge case) . Now we introduce 
another new bracket 

ai— r+l f3i—r+l 

(7.18) 

Then substituting ([7T6l) - (|7Tf|) into ([5TT4l) - ([5T5l) and using ([7T8l) . we obtain 



^-^ = {p.HpUys{q\pr)]„. (7.20) 
This new bracket governs the time evolution in the gauge case 

^ = {X,iI,.,.(9^pO}„- (7-21) 

The brackets (|7.5p and (|7.18p are not anticommutative and do not satisfy the 
Jacobi identity. Therefore, the standard quantization scheme is not applicable 
here directly. We expect that some more intricate further assumptions should 
make it possible to quantize consistently degenerate Lagrangian systems within the 
suggested approach (see, e.g., UHISlj). 

It is worthwhile to consider the limit case, when rj = 0, i.e. 

{q^.p) = (7.22) 

identically, which can mean that Ba {q'^TPi) = 0, so the Lagrangian can be inde- 
pendent of the degenerate velocities v". It follows from ()5.16p that 

DaHpnys {q\p.) = ^^^''^^Jf^'P') ^ 0, (7.23) 

which leads to the "independence" statement: the "physical" Hamiltonian 
Hphys {q^:Pi) does not depend on the degenerate coordinates q" iff the Lagrangian 
does not depend on the velocities v"". In this case, the bracket (|7.18p coincides with 
the Poisson bracket in the reduced ("physical" space) { , }gauge — { ' }phys- 
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Example 7.1 (Christ-Lee model [H]). The Lagrangian of SU (2) Yang-Mills theory 
in + 1 dimensions is (in our notation) 

L{Xi,ya,Vi) ^ ^ X! I ~ X! ^ija^jVa] ^ U {x'^) , (7.24) 
i=l,2,3 \ j,Q = l,2,3 / 

where i,a — 1,2,3, x^ — '^^xf, Vi — Xi and Sijk is the Levi-Civita symbol. Be- 



cause (|7.24l) is independent of degenerate velocities y^, all Ba (q^^Pi) 0, and 

therefore Fa/s {q^,Pi) Q, we have the limit gauge case of the above classifi- 

cation. The corresponding Clairaut equation (13. 3p for H = i/*-" {xi,ya,Pi,Pa) has 
the form 



1=1,2,3 Q=l,2,3 1=1,2,3 \ j,a=l,2,3 / 

(7.25) 

We show manifestly, how to obtain the envelope solution for regular variables and 
general solution for degenerate variables. Its general solution is 

Hgen = ^ Pid + ^ paCa " 2 X! \ ^ X! ^ijaXjya | + U {x"^) , 
i=l,2,3 a=l,2,3 j=l,2,3 \ j",a=l,2,3 J 

(7.26) 

where Ca are arbitrary functions of coordinateso We differentiate f|7.26p by 



- X! ^ija^jVa , (7.27) 
j,Q = l,2,3 / 

(7.28) 

and observe that only the first relation (I7.27P can be resolved with respect to 
and therefore can lead to the envelope solution, while other Cq cannot be resolved, 
and therefore we consider only general solution of the Clairaut equ ation . So we can 

exclude half of the constants using (|7.27l) (with the substitution pi pi) and get 
the mixed solution (|4.4p to the Clairaut equation (|7.25l) as 



■f^mL (a;»,ya,P»,PQ,Ca) = I ^ pf + ^ eijap.iXjya+ ^ PaCa + U{x^). 

(7.29) 



2 

1=1,2,3 i,j,a=l,2,3 Q=l,2,3 



Using (|5.10l) . we obtain the "physical" Hamiltonian 



{x^,ya,Pi) ^ \ ^ Pi + X! ^ijaPiXjya + U {x^) . (7.30) 



i=l,2,3 i,i,Q=l,2,3 



From the other side, the Hessian of (|7.24l) has rank 3, and we choose Xi, Vi and 
ya to be regular and degenerate variables respectively. The degenerate velocities 
Va = 2/a cannot be defined from (|5.16p at all, they are arbitrary, and the first 



■^Recall that are passive variables under the Legendre transform. 
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integrals (|5.16|) . (|7.23p of the system (|5.14|) - ()5.15p become (also in accordance to 
the independence statement) 



dHpiiys {Xi , ya 5 Pi 



SijaPiXj = 0. (7.31) 



The preservation in time (|7.2ip of (I7.3ip is fulfilled identically due to the antisym- 
metry properties of the Levi-Civita symbols. It is clear that only 2 equations from 
3 of (I7.3ip are independent, so we choose piX2 = P2X1, Pix^ = psXi and insert in 
([7:301) to get 

Hphys = lpl^ + U{x^). (7.32) 



The transformation p = piVx'^/'xi, x = V x^ gives the well-known result [11 ] [27 ] 

Hphys ^Ip^ + U (5) . (7.33) 

8. Singular Lagrangian systems and many-time dynamics 

The n-time classical dynamics and its connection with constrained systems were 
studied in [35l[36] as a generalization of some relativistic two-particle models [H] . 
Under the Clairaut-type approach, we now treat the degenerate coordinates as 
parameters analogous to n — r time variables (with n — r corresponding "Hamilto- 
nians" —Ba [q^,Pi), see (15.12^ 1. Indeed, let us introduce n — r + 1 "times" t** and 
the corresponding "many-time Hamiltonians" (t*^, q'',pij, fi — 0, .. .n ~ r hy 

t^tO, Hphys {q\q",P^)^J^o{q\t'',P^), (8.1) 
q^t^^a^ ^Ba,{q\q'',p,)^M'„{q\t'',Pi), a^l,...,n-r. (8.2) 

Then the equations (|5.14p - f)5.15p can be presented in the differential form 

n — r 

rf9' = E{9^^M}*^ (8-3) 

71— r 

dp, = ^fe,^^}dt^, (8.4) 
The system of equations for degenerate velocities (j5.16p becomes 

n — r 

J2 %udt'' = 0, (8.5) 

IJ,=0 



where 



(8.6) 



dt- ^t^^ ' ^' 
If we introduce a 1-form uj — Pidq^ — J^dt^, then it follows from (18. 5p that 



n — r n — r 

duJ^-J2J2 ^M''*^ ^ = 0, (8.7) 



2 

^=0 u=0 
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which agrees with the action principle in the form S = J uj. The corresponding set 
of the Hamilton- Jacobi equations for 5 = 5 (t, q^) i — > (t^, q*) is 

Therefore we come to the conclusion that any singular Lagrangian theory (in the 
Clairaut-type formulation |19j ) is equivalent to the many-time classical dynamics 
[T71I35]: the equations of motion are (|8.3l) - (|8.4p and the integrability condition is 
= 0. 

9. Conclusions 

To conclude, we have described Hamiltonian evolution of singular Lagrangian 
systems using n — r + 1 functions Hphys {q^,Pi) and Ba {^q"^,Pi) of dynamical vari- 
ables. This is done by means of the generalized Legendre-Clairaut transform, that 
is by solving the corresponding multidimensional Clairaut equation without intro- 
ducing Lagrange multipliers. All variables are set as regular or degenerate accord- 
ing to the rank of the Hessian matrix of Lagrangian. We consider the restricted 
"physical" phase space formed by the regular coordinates and momenta pi only, 
while degenerate coordinates q°' play a role of parameters. There are two reasons 
why degenerate momenta pa corresponding to need not be considered in the 
Clairaut-type formulation: 

1) the mathematical reason: there is no possibility to find the degenerate ve- 
locities ii", as can be done for the regular velocities in (|4.2I) . and the "pre- 
Hamiltonian" (12.21) has no extremum in degenerate directions; 

2) the physical reason: momentum is a "measure of movement", but in "de- 
generate" directions there is no dynamics, hence — no reason to introduce the 
corresponding "physical" momenta at all. 

Note that some possibilities to avoid constraints were considered in a different 
context in [T31I11] and for special forms of the Lagrangian irQ 25, . 

The Hamiltonian form of the equations of motion (|7.3I) - (|7.4I) is achieved by intro- 
ducing new brackets (|7.5p and (|7.18|) which are responsible for the time evolution. 
However, they are not anticommutative and do not satisfy the Jacobi identity. 
Therefore quantization of such brackets requires non-Lie algebra methods (see, 

e.g., mm)- 

In the "nonphysical" coordinate subspace, we can formulate some kind of a non- 
abelian gauge theory, such that nonabelianity appears due to the Poisson bracket 
in the physical phase space. 

Finally, we show that, in general, a singular Lagrangian system in the Clairaut- 
type formulation [Tnidn] is equivalent to the many-time classical dynamics. 
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Appendix A. Multidimensional Clairaut equation 

The multidimensional Clairaut equation for a function y = y{xi) of n variables 
Xi is 3, 32^ 

n 

y = (A.l) 

where prime denotes a partial differentiation by subscript and / is a smooth func- 
tion of n arguments. To find and classify solutions of (jA.ip . we need to find first 
derivatives y'^. in some way, and then substitute them back to (jA.ip . We differen- 
tiate the Glairaut equation (jA.ip by Xj and obtain n equations 

n 

1=1 

The classification follows from the ways the factors in (|A.2p can be set to zero. 
Here, for our physical applications, it is sufficient to suppose that ranks of Hessians 
of y and / are equal 

ranky"^^^. = rank/^'^^^^^ = r. (A.3) 

This means that in each equation from (|A.2p cither the first or the second multiplier 
is zero, but it is not necessary to vanish both of them. The first multiplier can be 
set to zero without any additional assumptions. So we have 

1) The general solution. It is defined by the condition 

y'L., = 0. (A.4) 

After one integration we find y^. ~ Ci and substitution them to (jA.ip and obtain 

n 

Vgen = ^ XjCj - f {d) , (A.5) 

where Ci are n constants. 

All second multipliers in (|A.2p can be zero for i = 1, ... ,71, but this will give 
a solution, if we can resolve them under y'^.. It may be possible, if the rank of 
Hessians / is full, i.e. r — n. In this case we obtain 

2) The envelope solution. It is defined by 

X^=fy'^^■ (A.6) 

We resolve (|A.6|) under derivatives as y^. = Q {xj) and get 

n 

Venv = ^XiC\ (Xj) ~ f{Ci (Xj)), (A. 7) 

i=l 

where Ci (xj) are n smooth functions of n arguments. 
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In the intermediate case, we can use the envelope solution (jA.7|) for first s vari- 
ables, while the general solution (jA.SP for other n — s variables, and obtain 
3) The s-mixed solution, as follows 

s n 

y^L = ^^^3 (^'j) + ^i^i ~ -^(^1 (^j) ' ■ • ■ ' C's {xj) , . . . , c„). (A.8) 

If the rank r of Hessians / is not full and a nonsingular minor of the rank r is 
in upper left corner, then we can resolve first r relations (jA.6p only, and so s < r. 
In our physical applications we use the limited case s = r. 

Example A.l. Let f (zi) = + z| + ^3, then the Clairaut equation for y = 

yixi,X2X3) is 

y = Xly'x^ + X2y'^^ + x^y'^,^ - {y'^y - {y'^y - y'^^, (A.9) 

and we have ti = 3 and r = 2. The general solution can be found from (jA.4[) by 
one integration and using (jA.sp 

ygen = Cl (xi - Ci) + C2 (X2 - C2) + C3 (X3 - 1) , (A. 10) 

where Ci are constants. 

Since r = 2, we can resolve only 2 relations from (jA.6l) by y'x^ = ^ ^ y'x2 — ^ ■ 
there is no envelope solution (for all variables), but we have several mixed solutions 
corresponding to s = 1, 2: 

M^ I -r + C2(2:2-C2)+C3(x3-l), 

2/iL=<! 4 ^2 (A.ll) 

Cl (a;i - Cl) + + C3 (2:3 - 1) , 

= f + 1+^3 (-3-1). (A.12) 

The case / [zi) — zf + Z2 can be obtained from the above formulas by putting 
a;3 = C3 = 0, while y^]^ becomes the envelope solution y^nv = ^ + 



Appendix B. Correspondence with the Dirac approach 

The connection of the Clairaut-type formulation with the Dirac approach is 
made by interpretation of the parameters pa entering to the general solution of 
the Clairaut equation as the "physical" degenerate momenta Pa using for them the 
same expression through the Lagrangian as (|4.2p 



Pa=Pa^^y^. (B.l) 
Then we obtain the primary Dirac constraints (in the resolved form and our notation 

(El) 

*a(g'^,PA) =0, a = r + l,...n, (B.2) 

which are defined now on the fuU phase space T*M. Using ([STO]) and ([STT]) . we can 
arrive at the complete Hamiltonian of the first-order formulation ^24] (corresponding 
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to the total Dirac Hamiltonian [T5] ) 

n 

= Hphys{q^,P^)+ «"*a(g^,PA), (B.3) 

a—r+1 

which is equal to the mixed Hamilton-Clairaut function (|4.4p with the substitution 
(jB.ll) and use of (|B.2[) . Then the Hamilton-Clairaut system of equations (|5.14p - 
(|5.15p coincides with the Hamilton system in the first-order formulation [24 

q^ = {q^,HT}^,^^l, PA^{pA,HT}f,ll^ = 0, (B.4) 

and (|5.16D gives the second stage equations of the Dirac approach 

71 

{^a,HT}f^U^{^c.,Hp,,y,}^^^^+ {*",*/3}/„,j«^=0, (B.5) 

I3=r+1 

where 

is the (full) Poisson bracket on the whole phase space T*M. Note that 

={*„,*/3}^„j,, (B.7) 

DaHphys {q^,Pi) ~ {^a,Hphys} f^ii ■ (B-8) 

Our cases 2) and 1) of Section [7] work as counterparts of the first and the second 
class constraints in the Dirac classification [15], respectively. The limit case with 
zero "g"-field strength" = (fT^ (see ((B7)) ) corresponds to the 

Abehan constraints ^51157] . 
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